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Parameter-free effective field theory calculation 
for the solar proton-fusion and hep processes 
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Spurred by the recent complete determination of the weak currents in two-nucleon systems up to 
0(Q 3 ) in heavy-baryon chiral perturbation theory, we carry out a parameter- free calculation of the 
threshold S-factors for the solar pp (proton-fusion) and hep processes in an effective field theory that 
combines the merits of the standard nuclear physics method and systematic chiral expansion. The 
power of the EFT adopted here is that one can correlate in a unified formalism the weak-current 
matrix elements of two-, three- and four-nucleon systems. Using the tritium /3-decay rate as an 
input to fix the only unknown parameter in the theory, we can evaluate the threshold S factors 
with drastically improved precision; the results are S pp (0) = 3.94 x (1 ± 0.004) x 10~ 25 MeV-b and 
Shep(0) = (8.6 ± 1.3) x 10 -20 keV-b. The dependence of the calculated S-factors on the momentum 
cutoff parameter A has been examined for a physically reasonable range of A. This dependence 
is found to be extremely small for the pp process, and to be within acceptable levels for the hep 
process, substantiating the consistency of our calculational scheme. 

PACS number: 12.39.Fe 24.85. +p 26.65. +t 



I. INTRODUCTION 

The standard approach to nuclear physics Q an- 
chored on wavcfunctions obtained from the Schrodingcr 
(or Lippman-Schwinger) equation with "realistic" phe- 
nomenological potentials has scored an impressive quan- 
titative success in describing systems with two or more 
nucleons, achieving in some cases accuracy that defies 
the existing experimental precision. We refer to this ap- 
proach as SNPA (standard nuclear physics approach). 
The advent of quantum chromodynamics (QCD) as the 
theory of strong interactions raises a logical question: 
What is the status of SNPA in the context of the fun- 
damental theory, QCD ? Put more bluntly, is SNPA (de- 
spite its great success) just a mo del- dependent approach 
unrelated to the fundamental theory ? In our view this is 
one of the most important issues in nuclear physics today. 
In this paper, we propose that SNPA can be properly 
identified as a legitimate component in the general ed- 
ifice of QCD. The next important question is: If SNPA 
is indeed a bona- fide element of QCD, how can we es- 
tablish an expansion scheme which includes SNPA as an 
approximation and which allows a systematic calculation 



of correction terms with error estimation ? This system- 
atic correction with an error estimation is not feasible 
with SNPA alone. 

Broadly speaking, formulating nuclear physics calcula- 
tions starting from effective field theories (EFTs) based 
on QCD calls for "double-step decimation" To start 
with, the "bare" EFT Lagrangian for strong interaction 
physics in the nonperturbative sector is defined at the 
chiral scale A x ~ 47r/„. ~ 1 GeV, with the parameters in 
the Lagrangian determined by suitably matching them 
to QCD at that scale ||. This Lagrangian cannot be 
directly used for studying complex nuclei. In order to 
render it applicable to low-energy nuclear physics, it is 
desirable, if not indispensable, that a "decimation" be 
made from the chiral scale A x down to what might be 
called the "Fermi-sea scale" Aps ~ fep, where kp the 
nuclear Fermi momentum. An effective Lagrangian re- 
sulting from the decimation to this Fermi-sea scale is 
expected to contain parameters that reflect "intrinsic" 
density-dependence (as suggested in, e.g.,Ref. Q) in ad- 
dition to the usual dense loop effects. The next stage 
of decimation consists of integrating out the effective de- 
grees of freedom and modes from the scale Aps to Aq ~ 



MeV. This stage corresponds to doing shell-model type 
calculations in finite nuclei @] and to formulating Fermi- 
liquid theory in nuclear matter || . The physics of heavy 
nuclei or nuclear matter will involve both decimations 
but, for light nuclei, one can bypass the second decima- 
tion and work directly with the chiral Lagrangian. 

The aim of this article is to describe a formalism that 
combines the high accuracy of SNPA and the power of 
EFT to make totally parameter-free predictions for elec- 
troweak transitions in light nuclei. To be concrete, we 
shall consider the following two solar nuclear fusion pro- 
cesses: 

pp : p + p^d + e + + v e , (1) 
hep : p + 3 He -> 4 He + e+ + v e . (2) 

We stress that in our EFT approach these processes in- 
volving different numbers of nucleons can be treated on 
the same footing. A concise account of the present study 
was previously given in ]tJ for the pp process, and in ||] 
for the hep process. 

The reactions (|l|) and (Q) figure importantly in as- 
trophysics and particle physics; they have much bearing 
upon issues of great current interest such as, for exam- 
ple, the solar neutrino problem and non-standard physics 
in the neutrino sector. Since the thermal energy of the 
interior of the Sun is of the order of keV, and since no ex- 
perimental data is available for such low-energy regimes, 
one must rely on theory for determining the astrophys- 
ical S-f actors of the solar nuclear processes. Here we 
concentrate on the threshold S-factor, 5(0), for the reac- 
tions (|]) and (||) . The necessity of a very accurate esti- 
mate of the threshold S'-factor for the pp process, S PP (Q), 
comes from the fact that pp fusion essentially governs 
the solar burning rate and the vast majority of the solar 
neutrinos come from this reaction. Meanwhile, the hep 
process is important in a different context. The hep reac- 
tion can produce the highest-energy solar neutrinos with 
their spectrum extending beyond the maximum energy 
of the 8 B neutrinos. Therefore, even though the flux of 
the hep neutrinos is small, there can be, at some level, 
a significant distortion of the higher end of the 8 B neu- 
trino spectrum due to the hep neutrinos. This change 
can influence the interpretation of the results of a re- 
cent Super-Kamiokande experiment that have generated 
many controversies related to neutrino oscillations P,ffo| . 
To address these issues quantitatively, a reliable estimate 
of S/ lep (0) is indispensable. 

The primary amplitudes for both the pp and hep pro- 
cesses are of the Gamow- Teller (GT) type (A J = 1, no 
parity change). Since the single-particle GT operator is 
well known at low energy, a major theoretical task is 
the accurate estimation of the meson-exchange current 
(MEC) contributions. 

The nature of the specific challenge involved here can 
be elucidated in terms of the chiral filter picture. If the 



MEC in a given transition receives an unsuppressed con- 
tribution from a one-soft-pion exchange diagram, then 
we can take advantage of the fact that the soft-pion- 
exchange MEC is uniquely dictated by chiral symme- 
try Q and that there is a mechanism (called the chi- 
ral filter mechanism) that suppresses higher chiral-order 
terms pTp2] ] . We refer to a transition amplitude to which 
the chiral filter mechanism is applicable (not applica- 
ble) as a chiral-protected (chiral-unprotected) case. It 
is known that the space component of the vector current 
and the time component of the axial current are chiral- 
protected, whereas the time component of the vector cur- 
rent and the space component of the axial current are 
chiral-unprotected (see Appendix A). This implies among 
other things that the isovector Ml and axial-charge tran- 
sitions are chiral-protected JlJ,|l5| , but that the GT tran- 
sition is chiral-unprotected. This feature renders the esti- 
mation of the GT amplitude a more subtle problem; the 
physics behind it is that MEC here receives significant 
short-ranged contributions the strength of which cannot 
be determined by chiral symmetry alone. 

The difficulty becomes particularly pronounced for the 
hep process for the following reasons. First, the single- 
particle GT matrix element for the hep process is strongly 
suppressed due to the symmetries of the initial and fi- 
nal state wave functions. Secondly, as pointed out in 
Refs. fL§] (referred to as "CRSW91") and (referred 
to as "SWPC92"), the main two-body corrections to the 
"leading" one-body GT term tend to come with the oppo- 
site sign causing a large cancellation. A recent detailed 
SNPA calculation by Marcucci et al. |jl8| , hereafter re- 
ferred to as MSVKRB, has re-confirmed the substantial 
cancellation between the one-body and two-body terms 
for the hep GT transition. The two-body terms there- 
fore need to be calculated with great precision, which is 
a highly non-trivial task. Indeed, an accurate evaluation 
of the hep rate has been a long-standing challenge in nu- 
clear physics [^9|. The degree of this difficulty may be 
appreciated by noting that theoretical estimates of the 
hep S'-factor have varied by orders of magnitude in the 
literature. 

As mentioned, in obtaining accurate estimates of the 
GT transition amplitudes, it is imperative to have good 
theoretical control of short-distance physics. We expect 
that a "first-principle" approach based on effective field 
theory (EFT) will provide a valuable insight into this is- 
sue. We therefore adopt here the approach developed 
in Refs. [ ^2|pC| ], which purports to combine the highly 
sophisticated SNPA with an EFT based on chiral dy- 
namics of QCD. Our starting point is the observation 
that, to high accuracy, the leading-order single-particle 



A more modern and complete discussion on this observation 
has recently been given by Ananyan, Serot and Walecka [Q. 



operators in SNPA and EFT are identical, and that their 
matrix elements can be reliably estimated with the use 
of realistic SNPA wave functions for the initial and final 
nuclear states. Next, we note that in EFT the oper- 
ators representing two-body corrections^ to the leading- 
order one-body term can be controlled by systematic chi- 
ral expansion in heavy-baryon chiral perturbation theory 
(HBxPT) ^0|. Then, since the ratio of a two-body ma- 
trix element to the leading-order one-body matrix ele- 
ment can be evaluated with sufficient accuracy with the 
use of the realistic SNPA wave functions^], we are in a 
position to obtain a reliable estimate of the total (one- 
body -I- two-body) contribution. This approach takes full 
advantage of the extreme high accuracy of the wave func- 
tions achieved in SNPA while securing a good control of 
the transition operators via systematic chiral expansion. 
For convenience, we will refer to this method, which ex- 
ploits the powers of both SNPA and EFT, as "MEEFT" 
(short for more effective EFT). MEEFT - which is close 
in spirit to Weinberg's original scheme [pTJ based on the 
chiral expansion of "irreducible terms" - has been found 
to have an amazing predictive power for the n+p — > rf+7 
process and several other processes ^3|. An al- 

ternative approach, which however is in line with our 
reasoning, has been discussed by Ananyan, Serot and 
Walecka @. 

An early HBxPT study of the pp process was made in 
Ref. Q (hereafter referred to as PKMR98) by four of the 
authors. The calculation in PKMR98 was carried out up 
to next-to-next-to-next-to- leading order (N 3 LO) in chi- 
ral counting (see below). At N 3 LO, two-body meson- 
exchange currents (MEC) begin to contribute, and there 
appears one unknown parameter in the chiral Lagrangian 
contributing to the MEC. This unknown constant, called 
d R in Ref. pi| , represents the strength of a four-nucleon- 
axial-current contact interaction. In Ref. [ p4f , since no 
method was known to fix the value of d R , the ci^-term 
was simply ignored by invoking a qualitative argument 
that the short-range repulsive core would strongly sup- 
press its contribution. Due to uncertainties associated 
with this argument, Ref. was unable to corroborate 
or exclude the result of the latest SNPA calculation [^5) , 
$2B = 0.5 ~ 0.8 %, where 623 is the ratio of the con- 
tribution of the two-body MEC to that of the one-body 
current (see below). 

The situation can be greatly improved by using 
MEEFT. As first discussed in Refs. f?],|] and as will be ex- 



pounded here, the crucial point is that exactly the same 
combination of counter terms that defines the constant 
d R enters into the Gamow- Teller (GT) matrix elements 
that feature in pp fusion, tritium /3-decay, the hep pro- 
cess, /^-capture on a deuteron, and v-d scattering and 
that the short-range interaction involving the constant 
d R is expected to be "universal," that is, A-independent. 
Therefore, assuming that three- and four-body currents 
can be ignored (which we will justify a posteriori), if the 
value of d R can be fixed using one of the above pro- 
cesses, then we can make a totally parameter-free predic- 
tion for the GT matrix elements of the other processes. 
Indeed, the existence of accurate experimental data for 
the tritium /3-decay rate, T%, and the availability of ex- 
tremely well tested realistic wave functions for the A=3 
nuclear systems allow us to carry out this program. In 
the present work we determine the value of d R from F^ 
and perform parameter-free EFT-based calculations of 
Spp(0) and S hep (0). 

As described below, our scheme has a cutoff parame- 
ter A, which defines the energy/momentum cutoff scale 
of EFT below which reside the chosen explicit degrees 
of freedom^. Obviously, in order for our result to be 
physically acceptable, its cutoff dependence must be un- 
der control. In our scheme, for a given value of A in a 
physically reasonable range (to be discussed later), d R 
is determined to reproduce T^; thus d R is a function 
of A. According to the premise of EFT, even if d R it- 
self is A-depcndcnt, physical observables (in our case the 
S-factors) should be independent of A as required by 
renormalization-group invariance. We shall show that 
our results meet this requirement to a satisfactory de- 
gree. The robustness of our calculational results against 
changes in A allows us to make predictions on S pp (0) 
and Shepify with much higher precision than hitherto 
achieved. Thus we predict: S pp (0) = 3.94 x (1 ± 0.004) x 
1(T 25 MeV-b and S hep (0) = (8.6 ± 1.3) x 1CT 20 keV-b. 

The remainder of this article is organized as follows. 
In Section II we briefly explain our formalism; in par- 
ticular, we describe the relevant transition operators de- 
rived in HBxPT. Section III presents the calculation of 
S pp (0), while Section IV is concerned with the estimation 
of Shep(0)- Section V is devoted to discussion and conclu- 
sions. We have made the explanation of the formalism in 
the text as brief and focused as possible, relegating most 
technical details to the Appendices. 



2 The argument made here should apply generally to n-body 
currents (n > 2) but since the 2-body terms are dominant, 
we shall continue to restrict our discussion to the latter. 

3 This statement holds only for the finite-range part of two- 
body operators, with the zero-range part requiring a regular- 
ization to be specified below. 



4 The cutoff specifies not just the relevant degrees of freedom 
but also their momentum/energy content. This should be un- 
derstood in what follows although we do not always mention 
it. 



II. FORMALISM 

We sketch here the basic elements of our formalism. 
The explicit degrees of freedom taken into account in 
our scheme are the nucleon and the pion, with all other 
degrees of freedom (p- and w-mesons, A(1232), etc.) in- 
tegrated out. The HB^PT Lagrangian can be written 
as 



with the chiral order A defined as 



X=d+e+- -2 



(3) 



(4) 



where d, e and n are, respectively, the numbers of deriva- 
tives (the pion mass counted as one derivative), external 
fields and nucleon lines belonging to a vertex. Chiral 
symmetry requires A > 0. The leading-order Lagrangian 
is given by 

C =B[ivD + 2ig A S • A] B - ~ £ C A (BT A B) 2 



f*Tr(iAHA fi ) + lfTr(x + ) , 



(5) 



where B is the nucleon field in HBxPT; g A — 1.2670 ± 
0.0035 is the axial- vector coupling constant and 
/tt = 92.4 MeV is the pion decay constant. Furthermore 

D fl B = (d^ + T fl )B, 

i 



with 



exp i 



2U 



(6) 



(7) 



R„ = £ (V° + Al) and = £ (V£ - At) denote ex- 
ternal gauge fields, and x is proportional to the quark 
mass matrix. If we neglect the small isospin-symmetry 
breaking, then x = (in the absence of external scalar 
and pseudo-scalar fields). For convenience, we work in 
the reference frame in which the four-velocity t> M and the 
spin operator are 



(1, 0) and = 



(8) 



The NLO Lagrangian (the so-called "1 /to" term) in the 
one-nucleon sector is given in Ref. |E6j, while that in 



the two-nucleon sector is given in Ref. [|3Q- With four - 
fcrmion contact terms included, the Lagrangian takes the 
form 



d=B 



v^v v - g^ v 
2tojv 



+ 2c 4 + 



2tojv 



D^Dv + 4c 3 zA • iA 



[£", S v ] [iA v , iA v ] 



.1J-C6 ^ s ^ j+^ B _ udl bs ■ AB BB 



m N 



+2id 2 e abc euu\sv fl A u > a BS X T h B BS s t c B 



(9) 



where tun — 939 MeV is the nucleon mass, and 

f»u = Zid^Ly - d v L^ - i [L M L u ])£ 

H j (d»Rv - d„R P -i[R„, R„])t, (10) 

£oi23 = 1, and A M = ^-A° We have shown here only 
those terms which are directly relevant to our present 
study. The dimensionless low-energy-constants (LECs), 
c's and gTs, are defined as 

C3,4 = 7— c 3 ,4, di,a = ___ gA J9 ^1,2- (11) 



m N 



m N f; 



We now consider the chiral counting of the electroweak 
currents (see the Appendices for details). In the present 
scheme it is sufficient to focus on "irreducible graphs" 
in Weinberg's classification pl[ . Irreducible graphs are 
organized according the chiral index v given by 



2{A - C) + 2L + J2 > 



(12) 



where A is the number of nucleons involved in the pro- 
cess, C the number of disconnected parts, and L the num- 
ber of loops; Vi is the chiral index A, eq.([|), of the z-th 
vertex. One can show that a diagram characterized by 
eq.(|l^) involves an n^-body transition operator, where 
hb = A — C + 1. The physical amplitude is expanded 
with respect to v. As explained at length in the Ap- 
pendix, the leading-order one-body GT operator belongs 
to v=0. Compared with this operator, a Feynman dia- 
gram with a chiral index v is suppressed by a factor of 
(Q / h x Y , where Q is a typical three-momentum scale or 
the pion mass, and A x <~ 1 GeV is the chiral scale.M In 



Our definition of the pion field here is different from that 
used in Ref. Q ; we have changed the sign of the pion field. 
Furthermore, we employ here manifestly Lorentz-invariant 
and chiral-invariant interactions. 

6 For convenience, a chiral order corresponding to v is of- 
ten referred to as NTO; v—\ corresponds to NLO (next-to- 
leading order), u—2 to N 2 LO (next-to-next-to-leading order), 
and so on. 



our case it is important to take into account the kine- 
matic suppression of the time component of the nucleon 
four-momentum. We note 



v ■ pi ~ v ■ p[ ~ v ■ ki 



m N 



(13) 



where pf (p/*) denotes the initial (final) momentum of 
the l-th nucleon, and k^ = (p[ — pi) 11 . Therefore, each 
appearance of v ■ pi, v ■ p[ or v ■ ki carries two powers 
of Q instead of one, which implies that v increases by 
two units rather than one. It is also to be noted that, if 
we denote by g M = (goi q) the momentum transferred to 
the leptonic pair in eqs.(|]) (||), then q ~ |g|~ Q 2 /A x 
~ 0(Q 2 ) rather than 0{Q) as naive counting would sug- 
gest. These features turn out to simplify our calculation 
considerably. 

In this paper, as far as the main calculation is con- 
cerned, we shall limit ourselves to N 3 LO; for certain dis- 
cussions, however, we shall consider operators belonging 
to N 4 LO as well. 

We now describe the derivation of one-body (IB) and 
two-body (2B) current operators with due consideration 
of chiral counting. The current in momentum space is 
written as 

J"(g)=V r "(Q)+A"(g) = J dxe- tqx J^(x). (14) 

When necessity arises to distinguish the space and time 
components of the currents, we use the notations 



V" = (V°, V), A" = (A°,A). 



(15) 



For the clarity of presentation, we first give a sum- 
mary chart of the basic chiral counting characteristics 
of the relevant currents, and then provide more detailed 
explanations in the remainder of this section and in the 
Appendices. The chiral counting of the electroweak cur- 
rents is summarized in Table [|, where the non- vanishing 
contributions at q — are indicated^] 

We now discuss the entries of this table order by order: 

• v = : One-body A and V°: A gives the Gamow- 
Teller (GT) operator, while V° is responsible for 
the charge operator. 

• v = 1 : One-body A and V: A gives the axial- 
charge operator while V gives the Ml operator. 



7 For small but finite values of q 7^ 0, there are slight devi- 
ations from this chart; for instance, for all the four cases in 
the table, there arise one-body contributions at NLO and at 
higher orders. However, these deviations are not significant 
in our case, which involves very small values of q. 



TABLE I. Contributions from each type of current at 
q — 0. The entry of "— " indicates the absence of contri- 
bution. "1B-RC" stands for relativistic corrections to the 
one-body operators, and "2B-1L" for one-loop 2-body contri- 
butions including counter term contributions. 



J" 


LO 


NLO 


N 2 LO 


N 3 LO 


N 4 LO 


A 


IB 




1B-RC 


2B 


1B-RC, 2B-1L and 3B 


A 




IB 


2B 


1B-RC 


1B-RC, 2B-1L 


V 




IB 


2B 


1B-RC 


1B-RC, 2B-1L 


v° 


IB 






2B 


1B-RC, 2B-1L and 3B 



• v = 2 : Two-body tree current with Vi — vertices, 
namely, the soft-pion-exchange current. This is a 
leading correction to the one-body Ml and axial- 
charge operators carrying an odd orbital angular 
momentum. 

• v = 3 : Two-body tree currents with J^. V{ = 1, 
which correspond to the hard-pion current, consid- 
ered in CRSW91 @ and SWPC92 @. These are 
leading corrections to the GT and V° operators car- 
rying an even orbital angular momentum. 

• v = 4 : All the components of the electroweak 
current receive contributions of this order. They 
consist of two-body one-loop corrections as well 
as leading-order (tree) three-body corrections. 
Among the three-body currents, however, there 
are no six-fermion contact terms proportional to 
(NN) 3 , because there is no derivative at the ver- 
tex and hence no external field. 

It is noteworthy that the counting rule for V is the 
same as for A , and the counting rules for V° and A 
are the same. The behavior of V and A summarized in 
Table | represents the chiral filter mechanism jllj , and V 
and A are referred to as chiral- filter-protected currents. 
By contrast, V° and A belong to chiral-filter-unprotected 
currents. 

We now discuss the explicit expressions for the relevant 
currents. For the one-body (IB) currents, for both the 
vector and axial cases, one can simply carry over the 
expressions obtained in MSVKRB. Up to N 3 LO, the IB 
currents in coordinate representaion are given as 



V°(l) =T- e - lq - r < 

V(0=rfe-^« 



1 + iq ■ cti x p l 



2n v - 1 



Pi 

m N 



Pi 

2m% 



Am 2 N 



■ W 

+ 1- — q x f 1 



2m 



N 



2pv - 1 

-101 x Pl q 



Am 



N 



A (I) = -jUTfe-*-" 



gj Pl 

m N 



Pl 

2m\ 



A(l) = -g AT ~e-^ r ' [a l 



2(P; (Ti-Pi- <Tipf)+iq x p l 



4m 2 N 



(16) 



where \iy ~ 4.70 is the isovector anomalous magnetic 
moment of the nucleon. Eq.(|T6|) gives the isospin- 
lowering currents, 



t _ ja= 



(17) 



and t 1 = i(rf — i-rf ). The tildes in eq.(|l6|) imply that 
the currents are given in the coordinate space represen- 



(Vj - Vi) 



Vi act on 



tation, and p l = — iVi and Pj = 
the wave functions. 

We next discuss the two-body (2B) currents. The ex- 
pressions for the V 2 b and A 2B operators can be found in 
|^,^0|. The V 2 ° B operator does not appear up to the or- 
der under consideration. The derivation of the 2B axial 
current, A2B, in HB^PT is described in Appendix A. In 
momentum space, A 2 b reads 



12B — 



l in • 



A12 



9A 



P [cri - a 2) 



+4 c 3 fc k ■ (tj cri + t 2 <x 2 ) 



' J + - J T x fe X \<J , 



X fe] 



9A 



2di(r 1 cri + t 2 <x 2 ) 



■ d 2 TxCr > 



with p = ( Pl - p 2 )/2, Pl = (p, + p/)/2, r, = i(rf - 
* r f )j r x = ( T i x T 2) x — i{ T i x T 2) y , and similarly for cr x ; 
c's and d's are the LECs explained in PKMR98. The 
values of c's in Eq.(18) have been determined from 7r- 
7V data §9|: c 3 = -3.66 ± 0.08 and c 4 = 2.11 ± 0.08. 
The two constants, d\ and d 2 , remain to be fixed but it 
turns out (see Appendix C.2) that, thanks to Fermi-Dirac 
statistics, only one combination of them, 



d R = di 



„i 2 ~ 1 

2d 2 + -C3 + -C4 + - 

3 3 6 



(19) 



is relevant in the present context |] 

It should be noted that the two-body currents given in 
eqs. ( ffg ) are valid only up to a certain cutoff A. This im- 
plies that, when we go to coordinate space, the currents 
must be regulated. This is a key point in our approach. 
Specifically, in performing Fourier transformation to de- 
rive the r-space representation of a transition operator, 



we use the Gaussian regularization (see Appendix C). 
This is, to good accuracy, equivalent to replacing the 
delta and Yukawa functions with the corresponding reg- 
ulated functions, 



VIM 

2/2AM 



d 3 k 



Sl(k 2 )e ik - r . 



(2tt) 3 
d 



k 1 



1 919 - 
ml r drrdr y ° A[r) - 



where the cut-off function, S\(k 2 ), is defined as 



5a (fe 2 ) = exp 



k z 

2 A 2 



(20) 



(21) 



The resulting r-space expressions of the currents in the 
center-of-mass (cm.) frame are 

Vi2(r) = -^-r- r [* x ■ <r 2 + S 12 i&(r)] 



9 A 



q x 



OMr) + ( T x - -0 X ) &(r) 



(18) A 12 (r) 
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2m N f 2 
+S A (r)d R 
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2/2aW 



(22) 



where r = T\ — r 2 , S12 = 3<T] • r <r 2 • f — cri ■ er 2 , and 

e>* = t-q4, e> = r^cr , Tq = ff ■ Oq — ±e> , 

= ±, x , Tq = (ri t 2 )~ = (ri t 2 ) x - i(n t 2 ) 9 and 
<tq = (<t 1 0<t 2 ). We emphasize again that A12 in Eq.([2^) 
contains only one unknown LEC, d R , that needs to be 
fixed using an empirical input. As mentioned in Section 
I, we choose here to determine d R using the experimental 
value of r« . 



III. DETERMINATION OF d R 



A sign error made in the expression for d R in PKMR98 
has been corrected here. 



The cutoff parameter A characterizes the energy- 
momentum scale of our EFT. A reasonable range of A 



may be inferred as follows. According to the general tenet 
of %PT, A larger than A x ~ 47r/ w ~ has no physical 
meaning. Meanwhile, since the pion is an explicit degree 
of freedom in our scheme, A should be much larger than 
the pion mass to ascertain that genuine low-energy con- 
tributions are properly included. These considerations 
lead us to adopt A = 500-800 MeV as a natural range. 

In the present work we use as representative values A 
= 500, 600 and 800 MeV, and for each of these values of 
A we adjust d R to reproduce the experimental value of 
Tp. With the use of the value of d R so determined, we 
evaluate the pp and the hep amplitudes.^] 

To determine d R from Fg, we calculate Fg from 
the matrix elements of the current operators evaluated 
for accurate A=3 nuclear wave functions. We em- 
ploy here the wave functions obtained in Refs. @J30| 
using the correlated-hyperspherical-harmonics (CHH) 
method |3l] , p^| . It is obviously important to maintain 
consistency between the treatments of the A=2, 3 and 4 
systems. We shall use here the same Argonne vis (AV18) 
potential j}3| for all these nuclei. For the A > 3 systems 
we add the Urbana-IX (AV18/UIX) three- nucleon poten- 
tial Jm] . Furthermore, we apply the same regularization 
method to all the systems in order to control short-range 
physics in a consistent manner. 

The values of d R determined in this manner are: 

d R = 1.00 ± 0.07 for A = 500 MeV , 

d R = 1.78 ±0.08 for A = 600 MeV, (23) 

d R = 3.10 ± 0.10 for A = 800 MeV , 

where the errors correspond to the experimental uncer- 
tainty in Fg. Once d R has been determined, we are in a 
position to make a parameter-free calculation of the tran- 
sition amplitudes for pp and hep, which will be described 
in the next two sections. 



IV. THE pp PROCESS 

It is convenient to decompose the matrix element of 
the GT operator into one-body and two-body parts 



M = Mib +M 



2B 



(24) 



We discuss them separately. In PKMR98, an extensive 
analysis was made of the leading-order (LO) one-body 



9 The masses of the light-quark vector mesons (p and u) are 
less than 800 MeV; therefore, with the use of A = 800 MeV, 
an accurate description of certain observables might require 
the explicit presence of the vector mesons, even though the 
probe energy /momentum in question is much smaller than 
800 MeV. Thus the results with the cutoff A = 800 MeV 
should be taken with some caution. 



matrix element M.^ N , with a focus on the connection 



between EFT and the effective range expansion. The 
results obtained with the AV18 potential |33[ were 



■M?b N = (IT 0.02 % T 0.07 % =p 0.02 %) 
x 4.859 fin , 



(25) 



where the errors are due to uncertainties in the scatter- 
ing length and effective ranges. The "full" one-body con- 
tribution in PKMR98 includes the vacuum-polarization 
(VP) and two-photon-exchange (C2) contributions. In 
PKMR98, howev er, t he one-body current due to the 
1/rn^, term in eq.( |A.3| ) was ignored. Although this term 
is required for formal consistency, its numerical value 



turns out to be quite small, M 



l/mi 



-0.006 fin. In 



terms of the A pp defined in Ref . ]35"| F] we have 



A , 



|a c | 2 7 3 



A%M 2 1B = 6.91 



(26) 



for the central value, where a is the pp 1 So scattering 
length, and 7 and As are the wave number and S-wave 
normalization constant pertinent to the deuteron, respec- 
tively. This should be compared with 6.93 obtained in 
Ref. @. 

The properly regularized two-body GT matrix ele- 
ments for the pp process read 



M 



2B 



dr 



{ 



m 2 ( 1\ 

-Z- I c 3 + 2c 4 + - j yJ A (r) u d (r) u pp (r) 



12r 



c 3 - c 4 



1 ) ylkir) w d{r) u P p(r) 



u d (r) - V2w d (r) u'(r) 



'd( r ) - ^2w'd( r )\ u p P ( r ) + —^-w d (r)u pp (r) 
d R df\r)u d (r)u pp (r)}, (27) 



where u d {r) and w d (r) are the S- and D-wave components 
of the deuteron wave function, and u pp (r) is the 1 So pp 
scattering wave (at zero relative energy) . The results are 
given for the three representative values of A in Table ||; 
for convenience, the values of d R given in Eq.(|2^) are also 
listed. The table indicates that, although the value of d R 
is sensitive to A, M2B is amazingly stable against the 
variation of A within the stated range. In view of this 



10 The subscript pp has been added here to avoid confusion 
with the cutoff parameter A. 



high stability, we believe that we are on the conservative 
side in adopting the estimate M2B — (0.039 ~ 0.044) fm. 
Since the leading single-particle term is independent of A, 
the total amplitude M. = M.\b + M-2B is A- independent 
to the same degree as M. 2 b- The A- independence of the 
physical quantity which is in conformity with the 
tenet of EFT, is a crucial feature of the result in our 
present study. The relative strength of the two-body con- 
tribution as compared with the one-body contribution is 

5 2B = 4^ = (0.86 ± 0.05) %. (28) 

MlB 

We remark that the central value of S2B here is consid- 
erably smaller than the corresponding value, 62B = 4 %, 
in PKMR98. Furthermore, the uncertainty of ±0.05 % 
in Eq.(p8|) is drastically smaller than the corresponding 
figure, ±4 %, in PKMR98. 

We now turn to the threshold S factor, S pp (0). Adopt- 
ing the value G v = (1.14939 ± 0.00065) x 10~ 5 GeV~ 2 
p6fl , we obtain 

S PP (Q) = 3.94 x ( I±*"lV f ( Ok) 2 

ppy ' \ 1.01 / V 1.2670/ I 6.91 j 

= 3.94 x (1± 0.0015 ±0.0010 ±e) (29) 

in units of 10~ 25 MeV-b. Here the first error is due to 
uncertainties in the input parameters in the one-body 
part, while the second error represents the uncertainties 
in the two-body part; e(~ 0.001) denotes possible un- 
certainties due to higher chiral order contributions (see 
below). To make a formally rigorous assessment of e, we 
must evaluate loop corrections and higher-order counter 
terms. Although an N 4 LO calculation would not involve 
any new unknown parameters, it is a non-trivial task. 
Furthermore, loop corrections necessitate a more elabo- 
rate regularization scheme since the naive cutoff regular- 
ization used here violates chiral symmetry at loop orders. 
(This difficulty, however, is not insurmountable.) These 
formal problems set aside, it seems reasonable to assess e 
as follows. We first recall that both tritium /3-decay and 
solar pp fusion are dominated by the one-body GT ma- 
trix elements, the evaluation of which is extremely well 
controlled from the SNPA as well as EFT points of view. 
Therefore, the precision of our calculation is governed by 
the reliability of estimation of small corrections to the 



dominant one-body GT contribution. Now, we have seen 
that the results of the present N 3 LO calculation nicely 
fit into the picture expected from the general premise of 
EFT: (i) the N 3 LO contributions are indeed much smaller 
than the leading order term; (ii) the physical transition 
amplitude M does not depend on the cutoff parameter. 
Although these features do not constitute a formal proof 
of the convergence of the chiral expansion used here, it 
is extremely unlikely that higher order contributions be 
so large as to completely upset the physically reasonable 
behavior observed in the N 3 LO calculation. It should 
therefore be safe to assign to £ an uncertainty compara- 
ble to the error estimate for the two-body part in Eq. (^9|) ; 
viz., e w 0.1 %. In this connection we remark that an 
axial three-body MEC contribution to the 3 H GT ma- 
trix element was calculated explicitly in SNPA jt8) and 
found to be negligible relative to the leading two-body 
mechanisms. This feature is consistent with the above 
argument since, in the context of EFT, the three-body 
MEC represents a higher-order effect subsumed in "e" in 
Eq.©. 

Apart from the notable numerical differences between 
the present work and PKMR98, it is worth noting that 
short-range physics is much better controlled in MEEFT. 
In the conventional treatment of MEC, one derives the 
coordinate space representation of a MEC operator by 
applying ordinary Fourier transformation (with no re- 
striction on the range of the momentum variable) to 
the amplitude obtained in momentum space; this cor- 
responds to setting A = 00 in Eq.@. In PKMR98, 
where this familiar method is adopted, the d R term ap- 
pears in the zero-range form, d R 8(r). PKMR98 chose 
to introduce short-range repulsive correlation with hard- 
core radius rc and eliminate the d R S(r) term by hand. 
The remaining finite-range terms were evaluated as func- 
tions of rc- M2B calculated this way exhibited substan- 
tial rc-dependence, indicating that short-range physics 
was not well controlled. Inclusion of the d R term, with 
its strength renormalized as described here, properly 
takes into account the short-range physics inherited from 
the integrated out degrees of freedom above the cutoff, 
thereby drastically reducing the undesirable (or unphys- 
ical) sensitivity to short-distance physics. 

V. THE hep PROCESS 



TABLE II. The strength d R of the contact term and the 
two-body GT matrix element, M2B, for the pp process calcu- 
lated for representative values of A. 



A (MeV) 


d R 


M2B (fm) 


500 


1.00 ±0.07 


0.076 - 0.035 d R ~ 0.041 ± 0.002 


600 


1.78 ±0.08 


0.097 - 0.031 d R ~ 0.042 ± 0.002 


800 


3.90 ±0.10 


0.129 - 0.022 d R ~ 0.042 ± 0.002 



In the notation of MSVKRB, the GT-amplitudes are 
given in terms of the reduced matrix elements, L\(q;A) 
and E\{q;A). Since these matrix elements are related 
to each other as E\{q;A) ~ y2Li(q;A), with the ex- 
act equality holding at q=0, we consider here only one 
of them, L\{q;A). For the three exemplary values of 
A, Table III gives the corresponding values of L\{q;A) 
at q = \q =19.2 MeV and zero cm. energy; for conve- 



nience, the values of d R in Eg. (p3|) are also listed. We 
see from the table that the variation of the two-body GT 
amplitude (row labelled "2B-total") is only ~10 % for the 
range of A under study. Note that the A-dependence in 
the total GT amplitude is made more pronounced by the 
drastic cancellation between the one-body and two-body 
terms, but this amplified A-dependence still lies within 
acceptable levels. 

TABLE III. Values of d R and Li(q;A) (in fm 3/2 ) for the 
hep process calculated as functions of the cutoff A. The in- 
dividual contributions from the one-body (IB) and two-body 
(2B) operators are also listed. 



A (MeV) 


500 


600 


800 


d H 


1.00 ±0.07 


1.78 ±0.08 


3.90 ±0.10 


Li{q;A) 


-0.032 


-0.029 


-0.022 


IB 


-0.081 


-0.081 


-0.081 


2B (without d R ) 


0.093 


0.122 


0.166 


2B (oc d R ) 


-0.044 


-0.070 


-0.107 


2B-total 


0.049 


0.052 


0.059 



TABLE IV. Contributions to S hep (0) (in 10~ 20 keV-b) 
from individual initial channels calculated as functions of A. 
The last column gives the results obtained in MSVKRB. 



A (MeV) 


500 


600 


800 


MSVKRB 




0.02 


0.02 


0.02 


0.02 


3 Si 


7.00 


6.37 


4.30 


6.38 


3 Po 


0.67 


0.66 


0.66 


0.82 


x Px 


0.85 


0.88 


0.91 


1.00 


3 A 


0.34 


0.34 


0.34 


0.30 


3 P 2 


1.06 


1.06 


1.06 


0.97 


Total 


9.95 


9.37 


7.32 


9.64 



Table |v| shows the contribution to the S-factor, at 
zero cm. energy, from each initial channel. For com- 
parison we also give the results of MSVKRB for the 
AV18/ULX interaction. It is noteworthy that for all the 
channels other than 3 Si, the A-dependence is very small 
(< 2 %). The 3 Si channel is the most sensitive to short- 
distance physics because the extraordinary suppression of 
the one-body GT contribution makes more pronounced 
the chiral-non-protected nature of the GT transition. In 
fact, the sensitivity of the 3 Si channel to short-distance 
physics would be larger if the contribution of the A term, 
which is rather sizable here despite its generic 1/m sup- 
pression, were omitted. It is therefore reassuring that 
the chiral-filter mechanism allows reliable estimation of 
the A term in this channel as well (besides the P-wave 
channels), see (jl8). 



Summarizing the results given in Table IV, we arrive 



at a prediction for the hep S-factoij^]: 

S hep {0) = (8.6 ± 1.3) x KT 20 keV-b , 



(30) 



where the "error" spans the range of the A-dependcncc 
for A=500-800 MeV. This result should be compared 
to that obtained by MSVKRB @, S hep (0) = 9.64 x 
KT 20 keV-bQ 

The latest analysis of the Super-Kamiokande data ]3^ j 
gives an upper limit of the solar hep neutrino flux, 
<Miep) SK < 40xl0 3 cm s . The standard solar model 
@ using the hep S-factor of MSVKRB |J| predicts 
$(hep) SSM = 9.4 x 10 3 cm _2 s -1 . The use of the cen- 
tral value of our estimate, Eq.(|30|), of the hep S-factor 
would slightly lower &(hep) SSM but with the upper limit 
compatible with <fr(hep) SSM in Ref. J38fl. A significantly 
improved estimate of Sfc, e p(0) in Eq.(|30|) is expected to 
be useful for further discussion of the solar hep problem. 

To reduce the uncertainty in Eq.(|3(]), we need to re- 
duce the A-dependence in the two-body GT term. Ac- 
cording to a general tenet of EFT, the cutoff dependence 
should diminish as higher order terms get included. In 



fact, the somewhat rapid variation seen in d R (Table III) 
and in the 3 Si contribution to S/j ep (0) (Table IV) as A 



approaches 800 MeV may be an indication that there is 
need for the explicit presence of the vector-mesons (p and 
uj) with mass my < A. This possible insufficiency could 
be remedied to a certain extent by going to higher orders. 
A preliminary study |]39|l indicates that it is indeed possi- 
ble to reduce the A-dependence significantly by including 
N 4 LO corrections. We expect that the higher order cor- 
rection would make the result for A = 800 MeV closer to 
those for A = 500, 600 MeV, bringing the MEEFT results 
closer to what was obtained in MSVKRB. This possibil- 
ity is taken into account in the error estimate given in 
Eq.©. 



VI. DISCUSSION 

It is worth emphasizing that the above MEEFT pre- 
diction for $2B for the pp process is in line with the latest 
SNPA results obtained in Ref. [0 (and mentioned ear- 
lier). There too, the short range behavior of the axial 
MEC was constrained by reproducing ri. The inher- 
ent model dependence of such a procedure within the 
SNPA context was shown to be very weak simply because 
at small inter-particle separations, where MEC contri- 
butions are largest, the pair wave functions in different 



11 See below for a possible caveat on the given error estimate. 

12 The earlier studies ]l6| , |l7| were based on less accurate 
variational wave functions than used here and in MSVKRB; 
furthermore they did not include P-wave capture contribu- 
tions, which account for w 40 % of the total S'-factor. 



nuclei are similar in shape and differ only by a scale fac- 
tor E3]. As a consequence, the ratios of GT and pp- 
capture matrix elements of different two-body current 
terms are nearly the same, and therefore a knowledge 
of their sum in the GT matrix element is sufficient to 
predict their sum in the pp-capture matrix element p5| . 

In order to better understand how the present scheme 
works, it is helpful to compare the hep reaction with 
the radiative np-capture. The polarization observables in 
ft + p — » d + 7 are known to be sensitive to the isoscalar 
Ml matrix element, MIS', and this amplitude has been 
extensively studied in EFT [ p2||4lj| . The similar features 
of the hep GT amplitude and MIS' are: (i) the leading 
one-body contribution is suppressed by the symmetries 
of the wave functions; (ii) there is no soft-pion exchange 
contribution; (iii) nonetheless, short-range physics can be 
reliably subsumed into a single contact term. In the ftp 
case the strength of this term can be determined from 
the deuteron magnetic moment (for a given value of A). 
The calculation in Ref. [^2| demonstrates that the A- 
dependence in the contact term and that of the remain- 
ing terms compensate each other so that the total MIS 
is stable against changes in A. This suggests that, if we 
go to higher orders, the coefficient of the contact term 
in question will be modified, with part of its strength 
shifted to higher order terms; however, the total physi- 
cal amplitude will remain essentially unchanged. These 
features are quite similar to what we have found here for 
the hep GT amplitude. 

Evaluating the matrix element of the leading-order 
one-body operator in EFT with the use of realistic nu- 
clear wave functions is analogous to fixing parameters in 
an EFT Lagrangian (at a given order) using empirical 
inputs j4^] ; the realistic wave functions in SNPA can be 
regarded as a theoretical input that fits certain sets of 
observables. In the present MEEFT scheme, we take the 
view that the same realistic wave functions also provide 
a framework for reliably calculating (finite-range) many- 
body corrections to the leading-order one-body matrix 
element. The short-ranged part inherited from the inte- 
grated out degrees of freedom is regulated by the d R term. 
This way of handling "short-range correlation" is analo- 
gous to Bogner et a/.'s derivation || of "V/ OUJ _fe" based 
on renormalization-group theory. While our approach 
here is, in certain cases, not in strict accordance with the 
systematic power-counting scheme of EFT proper, never- 
theless the severity of this potential shortcoming should 
vary from one case to another (see discussion in Ref. jl3| ) . 
For the pp and hep amplitudes under consideration, the 
degree of A-dependence exhibited by the numerical re- 
sults does suggest that deviations from rigorous power- 
counting cannot be too significative. Indeed, this type of 
"resilience" may also explain why the SNPA calculation 
in Ref. [fl8i gives a result very similar to the present one. 
It is true that the two-body terms in MSVKRB are not 
entirely in conformity with the chiral counting scheme we 



are using here; some terms corresponding to chiral orders 
higher than N 3 LO are included, while some other terms 
which are N 3 LO in EFT are missing (see Appendix B.3). 
Most importantly the (i^-term - that plays a crucial role 
here - is omitted in MKSVRB although heavy-meson ex- 
change graphs may account for some part of it. This for- 
mal problem, however, seems to be largely overcome by 
the fact that also in MSVKRB a parameter (the axial 
NA coupling strength) is adjusted to reproduce T^. 

Not unrelated to the above issue of power-counting 
is the question of consistency of embedding "realistic" 
wave functions obtained from "realistic" potentials that 
are fitted accurately to experiments into an EFT frame- 
work with the currents obtained to a given order of chiral 
perturbation theory. It is a well-known fact that poten- 
tials that fit experiments are not necessarily unique. The 
non-uniqueness resides however in the short-range part 
of the potential, with the long-range part primarily gov- 
erned by the pion exchange. Let us suppose that one 
can calculate potentials to a very high order in a consis- 
tent expansion (that is, consistent with symmetries etc.). 
The structure of the potential would depend on various 
aspects of the calculation. For instance, although they 
all may fit equally well various experimental data such as 
e.g., nucleon-nucleon scattering, different regularizations 
would lead to different potentials, the difference residing 
mainly in the short-range part. One might worry that 
this non-uniqueness would upset the basic premise of an 
EFT, rendering the predictions untrustworthy. 

Another intricate issue, which is also connected to 
short-range physics, is the off-shell ambiguity. This prob- 
lem should be absent in a formally consistent EFT. In 
MEEFT, however, we insert the current operators de- 
rived from irreducible diagrams up to a given chiral order 
between phenomenological (albeit realistic) wave func- 
tions. Since the inserted current involves off-shell par- 
ticles, there can in principle be terms other than those 
that have been included in our approach. While those 
additional terms that may be required to eliminate the 
off-shell dependence are expected to be of higher order 
than N 3 LO, this issue warrants a further examination. 

To answer the above question with full rigor, much 
more work is needed. However, partial and yet reason- 
ably satsifactory answers can be obtained from this work. 
For chiral-filter protected processes, we have presented 
a clear argument that the above-mentioned ambiguity 
does not matter at the level of accuracy in question. 



The results listed in Table IV for the P-wave capture (to 
which the chiral-protected time component of the axial 
current contributes) demonstrate this point. The ques- 
tion of short-distance ambiguity arises only for chiral- 
unprotected processes like the GT transition. As already 
explained, however, the d R regularization essentially re- 
moves this ambiguity. The point is that the physics of the 
degrees of freedom above the cutoff scale A gets lodged 



in the short-range d R term. In fixing this term as a func- 
tion of A via the experimental value of Tg, one is es- 
sentially incorporating the short-range correlations that 
render low-energy physics insensitive to short-distance 
physics. 

As for the off-shell problem, we note that for processes 
involving a long-wavelength external current - such as 
the solar pp and hep reactions - the off-shell ambiguity 
should be small, so long as one uses high-quality phe- 
nomenological wave functions that accurately describe 
processes without the external current. The wave func- 
tions used here describe with high accuracy a rich ensem- 
ble of data available for the systems in question; they de- 
scribe very well the three-nucleon scattering states, and 
furthermore, the n 3 He elastic scattering cross section as 
well as the coherent scattering length calculated with 
these wave functions are in excellent agreement with the 
experiments. What is involved here seems to be a generic 
feature. A similar stabilizing mechanism is at work when 
Bogner et al. g arrive at a unique effective force Vi ow —k 
by integrating out the high-energy/momentum compo- 
nents contained in various "realistic" potentials. Nuclear 
physics calculations done with this effective force Q] 
have much in common with the MEEFT calculation de- 
scribed here. Furthermore, we remark that different off- 
shell properties reflect different choices of the field vari- 
ables and that, for each choice, the LECs need to be 
readjusted. It is in principle possible to choose the field 
variables in such a manner that off-shell contributions be- 
come highly suppressed. We are essentially adopting this 
particular choice by using the forms of the transition op- 
erators described above and adjusting the corresponding 
LEC, d R , to reproduce Tt. 

A possible approach that is formally consistent with 
systematic power counting is the pionless EFT based 
on the power divergence subtraction (PDS) scheme pj| 
(for a recent review, see Ref. [^6|), which has been ap- 
plied to the pp fusion fl^ . Due to the fact that this 
scheme also involves one unknown low-energy constant, 
PDS has not so far led to a definite prediction on the pp 
fusion rate. The problem is that this approach cannot 
be readily extended to systems with A > 3, in partic- 
ular to electroweak transition amplitudes in these sys- 
tems. What is lacking presently is a method to correlate 
in a unified framework the observables in different nu- 
clei (different mass numbers). This limitation keeps one 
from exploiting the experimental data available for the 
A > 3 nuclei to fix unknown LEC. Apart from the basic 
problem of organizing chiral expansion for complex nuclei 
from "first-principles" , a plethora of parameters involved 
would present a major obstacle. (For recent efforts in this 
approach, see Ref. [j46||48| and references given therein.) 
This difficulty is expected to be particularly pronounced 
for the hep reaction. 

There has been a series of intensive studies by the 
Julich Group to extend EFT calculations in the Weinberg 



scheme to systems with three or more nucleons pj| . The 
relationship between this approach and the phenomeno- 
logical potential approach has been examined in great de- 
tail. This line of study, however, has been so far limited 
to nuclear observables that do not involve the electroweak 
currents. An extension of the formalism developed in 
Ref. |Q to electroweak transitions should be extremely 
useful. 
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APPENDIX A: GAMOW-TELLER OPERATORS 

The aim of this and subsequent Appendices is to pro- 
vide some technical details that have been left out in the 
main text. The readers who are not interested in the de- 
tails of our calculation can safely skip these Appendices 
without missing the essential points of our results. 

We decompose the axial current into n^-body opera- 
tors as 



A — /i 1B -f- A 2B -t /i 3B -f • 

= J2 A r + J2 A i£+ Y, A i™ 

/ j l / j hill / j LIILIL 



l<m /<m<n 



(A.l) 



where (I, m, n) are particle indices. The one-body oper- 
ator can be read from 



(N(p')\A^ a (x = 0)\N(p)) 
= -u{p') 



2m n 



-u(jp), (A.2) 



where u(p) is a four-component Dirac spinor of momen- 
tum p, and k ,J '=(p — p')** is the momentum carried by 
the lepton pair. Ga(A 2 ) and Gp(k 2 ) are the axial and in- 
duced pseudoscalar form factors, respectively. Note that 
k^=0(Q 2 /m N ), while p=0 (Q) and p '=0 (Q) . Thus the 
contribution from Gp term is of 0(Q 4 ). In getting the 
non-relativistic operators from the above relativistic form 
factors, we should also take into account the wave func- 
tion normalization. The resulting one-body operator up 
to 0(Q 3 ) then reads 
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This expression has appeared in eq.rtjl]). 

In the following subsections, we derive all the two-body 
GT operators up to N 4 LO and leading three-body GT 
operators . 



1. Two-body GT 

There are no two-body GT diagrams that involve only 
Vi = vertices, because the A 1 '°"kNN vertex is kinemat- 
ically suppressed, and there is no four-fermion contact 
contribution at LO (i/j = 0). As a consequence, the two- 
body GT operator starts at v = 3. The two-body GT 
operator at threshold (g M — + 0) was given up to N 3 LO 
in Ref. J50|. We extend here that analysis to include all 
contributions up to N 4 LO. To this end, it is useful to 



decompose Af m as 



A lm = A i m ( l7T ) + ^L( 27r )> 



(A.4) 



where A" m (l7r) represents the contributions of the one- 
pion pole part and A° m (27r) stands for the remaining 
short-ranged part. Generic diagrams for j4" m (l7r) and 
Af TO (27r) are shown in Fig. 1. 

We now list all the two-body GT operators belonging 
to v=i and v=A. 

• v = 3: 

This contribution comes from tree graphs (one-pion- 
exchange and contact) with a Vi = 1 vertex. The result- 
ing GT operators, given in Ref. |5C(], are of the form 
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(A.5) 



A^ 3 (2tt) = 



9A 
m N f 2 




FIG. 1. A one-pion pole diagram, (a), responsible for 
Af m (l7r), and a short-range contribution diagram, (6), re- 
sponsible for Af m (2n). The solid circles include counter-term 
insertions and (one-particle irreducible) loop corrections. The 
wiggly line stands for the external field (current) and the 
dashed line for the pion. One-loop corrections of the relevant 
orders for the pion propagator and the nNN vertex need to 
be included. 
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Although there are two unknown parameters, d\ and d 2 , 
it turns out that the Fermi-Dirac statistics effectively re- 
duces the number of unknowns to one. We will come 
back to this important point later. 

• v = 4: 

Tree graphs with v \ = 2 and one-loop graphs with 
J2i v i — enter at this order. Since there is no irNN 
vertex with i/j = 1, a v = 4 tree graph should have either 
A(NN) 2 or AttNN vertex with y t =2. We can, however, 
exclude either possibility. The absence of AttNN ver- 
tex at Vi=2 can be ascertained by consulting a complete 
list of terms that appear in the N 2 LO Lagrangian given 
in Ref. [|lj. As for the v i = 2 A(NN) 2 vertex for the 
two-nucleon sector, a complete list is not available yet. 
We therefore resort to parity selection rules. Our vertex 
must have one A M and one involving four nucleon 
fields. These operators should not be contracted with 
the four- velocity , for otherwise the actual chiral index 
would acquire an extra power of Q. We can easily show, 
however, that it is impossible to construct a parity-even 
Lorentz scalar with one A„, one and arbitrary num- 
bers of 5 M and ^ va ^ . Introducing an operator of the 
9fj,A v — d u A^ type instead of A M and does not help 
either. These observations lead us to conclude that no 
divergences occur in the relevant loops and, more impor- 
tantly, that no new parameters appear at i>=A. 

The one-pion-exchange contribution can be read off 
from 1-loop corrections to the vertex; the relevant 

diagrams are shown in Fig. |[ We note that only the first 
five diagrams (a)-(e) contribute to A. Using the expres- 
sions given in Ref. |lq| , where all the one- loop diagrams 
have been calculated, we find 
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FIG. 2. One loop diagrams that contribute to the A^nNN 
vertex. Only the first five diagrams (a) — (e) contribute to A. 
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where fc; = (I = 1, 2), and Di(k)'s are defined as 
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Here k = |fe|, z = 1 — z, and Afjfc = <J to 2 + zz k . 

The integrations over z can be done analytically, re- 
sulting in 
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with -f < 9 fc < f . 



Note that the one-pion-pole cont ribu tions can be ab- 
sorbed into A^ 3 (l7r) (given in Eq.( |A.5|) ) by renormaliz- 
ing £3 and £4, 
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For the two-pion contribution A1^ 4 (2it), the relevant 
one-loop graphs are shown in Fig. |^. Among the dia- 
grams in the figure, only the first four graphs, (a) — (d), 
can contribute to the GT operator; (e) is identically 
zero due to isospin symmetry, and the remaining graphs, 
(/) — (h), contribute only to A° . 
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FIG. 3. One loop diagrams for the A^ 4 (27r): the first 
four diagrams (a) — (d) contribute to the space part of the 
axial-vector (GT) and the remaining diagrams (/) — (h) to 
the axial-charge current. The graph (e) vanishes. 



As mentioned, the four diagrams (a) — (d) are all 
ultraviolet-finite. The first three graphs give 

A a 12 (2n : a) = --J^L.^ [3^(fc 2 ) + k 2 2 D 2 (k 2 )} 

+ (1^2), 

A a 12 {2n:b) = -^-j-^[a 1 D 1 {k 2 ) 

- (2k 2 (ri-k 2 -(Tikl)D 2 (k 2 )] 
+ (!«-> 2), 

A?2(27F : C) = " l^r/l^ + 2T * ] ^(k) 
+ (2fc 2 cr\ k 2 — <r\ k 2 )D^{k 2 ) 
+ (er 2 k\ - k 2 a 2 - k 2 )D 5 (k 2 )] 
+ (!.<->• 2) (A.13) 

while the fourth diagram (d) gives 

A a 12 (2n : d) = - fcj^ {CUT^, r^rM 
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FIG. 4. Diagrams for the three-body GT operator. All 
other diagrams are higher order than N 4 LO except for the 
crossed diagrams of (a) - crossed with respect to the particle 
indices. 



3. Comparison with SNPA exchange currents 



The summation here is taken over all possible combi- 
nations of spin-isospin operators (with no derivatives) 
that figure in the nucleon-nucleon interactions. Using 
a generic expression 



A 



Ca^A^A = X\ + <T\ ■ (J 2 X a + T\ ■ T 2 X T 



(Ti ■ er 2 fi ■ r 2 X aT (A. 15) 



where X\, X a , X T , X aT are constants that characterize 
the LO short-range nuclear forces, we can write 



A a 12 (2ir : d) = - 
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[(-aXi + Mf^XTfin+T^ra) 



327T/2 

+ {-2X aT ){fi x f 2 ) Q (<7i x er 2 ) 

+ {9X a - 3X T )(7f o-a + t^cti)] . (A.16) 

We demonstrate below that X's can all be absorbed into 
the parameters, d's. 

2. Three-body GT 

The three-body GT operators up to N 4 LO come from 
the two diagrams given in Fig. ^. They contain only 
vi = vertices, and their contributions read 
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The meson-exchange currents in SNPA |5^,|53|] are 
based on one-boson exchange diagrams involving those 
bosons which are responsible for the phenomenological 
nuclear forces in the context of one-boson-exchange mod- 
els. This framework does not have direct contact with 
chiral counting. We give here a detailed comparison be- 
tween the transition operators used in SNPA and those 
dictated by %PT. Among the most elaborate SNPA oper- 
ators are the ones used in CRSW91 these operators 
were derived by Towner [)53"| based on a phenomenologi- 
cal Lagrangian S] which satisfies CVC, PCAC and cur- 
rent algebra. We consider the SNPA operators used in 
CRSW91 as a representative. It will turn out that there 
are substantial differences between the SNPA and xPT 
operators in both the long-range and short-range parts. 

In CRSW91 the heavy particles, p and A, are treated 
as explicit degrees of freedom.^] To examine the roles of 
these heavy particles in the context of the present com- 
parison, we divide the two-body currents in CRSW91 
into two families: 



(A.17) 



A a = A a I +A a II (A.19) 
= [A a (A7r) + A a (np) + A a (TrS)} 

+ [A a (Ap)+A a ( P S)}, (A.20) 



13 The one-body operators are not sensitive to these addi- 
tional ingredients as long as the single-nucleon parameters are 
determined empirically. We note, however, that CRSW91 in- 
cludes only the leading-order GT operator, ^qac, without 
taking into account relativistic corrections. 



where the "S" stands for "seagull" . and A a u can be 
associated, respectively, with Af m (lTr) and Af TO (27r) in 
Eq.([A.4|). The expression for A a j is^] 
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(A.23) 



where mA ~ 1232 MeV, k ~ 6.6, and g p ~ 2.50 is the 
pNN coupling constant; A^ (A„) is a cutoff parameter 
characterizing the ttNN (pNN) coupling form factor. 
We have defined I\ and I 2 by 
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We note that I\ = 1 if we assume the Goldberger- 
Treiman relation, and / 2 = 1 if the KSRF relation 
holds .q| The above equation should be compared with 
A^ 3 (l7r) in Eq.QA.6;). A little exercise shows that, while 
the currents A(irA) and A(irp) can be related to certain 
currents in %PT, A(ttS) has no xPT counterpart to the 
order considered here. A possible explanation for the 
occurrence of this "extra" term in SNPA is that A(irS) 
arises as a "recoil" term associated with the use of the 
pseudoscalar coupling. A x?T analog of A(nS) would 
be a 1/mjv term, but this term should be substantially 
suppressed; hence a term like A(wS) should be absent in 
chirally invariant theory. Comparison of the coefficients 
of (T 1 XT 2 ) a p 1 , r^k 2 , (fixf 2 ) a cr 1 xfc2 and (fi x f 2 ) a <ri x q 
leads to the following correspondence between xPT (left- 
hand side) and CRSW91 (right-hand side): 



14 



The pseudovector ttNN coupling constant f^NN used in 
CRSW91 is related to the quantities used here as 



fnNN _ SttJVJV _ <?A 

m,r 2mjv 2/^ 



(A.21) 



For the last equality, the Goldberger-Treiman relation has 
been used. 

15 Since g p used in CRSW91 is half of the conventional one, 
the KSRF relation here reads 2(2g p ) 2 f% = m 2 , rather than 
2»p/^ = m%. 
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(1 + k). (A.28) 



The presence of the momentum-dependence in IZ's and 
the p-meson propagator prevents us from going beyond 
this correspondence. To proceed, however, we may con- 
sider the approximation 



1 « J 2 ^(fc 1 )72 w (fc 2 ) 
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We then find 
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It is informative to compare these values of c's with those 
obtained in a resonance-exchange saturation analysis by 
Bernard, Kaiser and Meissner (BKM) p6[ . We note that 
the two approaches give very different results for the A- 
resonance contributions. CRSW91 used the quark model 
value for the ratio, gANA/g^NA, the accuracy of which 
is rather difficult to assess. Meanwhile, the resonance- 
saturation calculation suffers from ambiguity related to 
the so-called off-shell parameter Z. Considering these un- 
certainties, it is perhaps not too surprising that BKM's 
estimate of the A contribution to 63, \c$ |=3.59, is 2.2 
times larger than the estimate in CRSW91. We also note 
that, while CRSW91 only includes the A and p- meson 
contributions, BKM's calculation contains the contribu- 
tions from scalar-meson and Roper exchanges as well. 
According to BMK, 



~scala 



gRoper = _ 3 _ _ Q Q6 = _ 4 



c^ + c p A + c4° por = 1.80 + 1.53 + 0.11 = 3.44. (A.33) 

Thus the contributions of the scalar-meson exchange are 
substantial. What is significant for our calculation is the 
fact that the coefficients c's can be extracted directly 
from the ttN scattering data |^,^l|. The most recent 
analysis |5lj] gives 

c 3 = (-5.58 ±0.08, -5.49 ±0.01, -5.82 ±0.08), 

c 4 = (3.26 ± 0.05, 3.29 ± 0.01, 3.30 ± 0.04). (A.34) 



These values are in reasonable agreement with those ob- 
tained in the resonance saturation approach. We should 



note, however, that the results in Eq.(A.34) belong to an 
N 4 LO calculation wherein A(1232) as well as other mas- 
sive degrees of freedom have been integrated out. The 
explicit inclusion of A(1232) would modify the values of 
c's, because the A contribution to the c's should now 
be excluded. We also should pay attention to a simi- 
lar modification of the LECs as we move from N 4 LO to 
N 3 LO. Although the difference between c's obtained in 
an N 3 LO calculation and those obtained in N 4 LO are of 
order of Q and hence can in principle be neglected in an 
N 3 LO calculation, it is more natural and safer to use in 
our present calculation the values obtained in an N 3 LO 
analysis f2q|, 
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c 4 



-4.96 ±0.23, 
3.40 ±0.09. 



(A.35) 



The precise value of Cq is unimportant in the present 
context, since it is suppressed by the kinematic factor 
\q\. In any event, the results of BMK and CRSW91 are 
close to each other, c§ ~ 5.83. 

We now discuss the short-ranged currents, A a u = 
A(pA) + A(pS). According to CRSW91, the dominant 
term in A a u is of the form 
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It should be noted, however, that this term belongs to 
N LO in chiral counting, and therefore its inclusion in 
CRSW91 constitutes a deviation from xPT. Although a 
particular N 5 LO term may give an appreciable contribu- 
tion (see below), there are many terms of the same or- 
der, including multi-loop diagrams, and in general there 
should be a substantial cancellation among these terms 
to make the net N 5 LO contribution small, as dictated by 
chiral symmetry. Thus, there are important differences 
between A a u of CRSW91 and A{2n) derived in X PT. 



the soft one-pion-exchange, which is of 0(Q 2 ). Loop con- 
tributions start at 0(Q 4 ), and hence the ratio of the loop 
contribution to the tree diagram two-body contribution 
is 0(Q 2 ). Finally, there is no three-body contribution up 
to 0(Q 4 ). 

The one-body axial-charge operator at threshold is 
given by 
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m N 
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(B.l) 



which is 0(Q 1 ). We note that there is no relativistic cor- 
rection of O(q) to the one-body axial-charge; this aspect 
is in sharp contrast to the GT operator. 

The two-body ^4° current appears at 0(Q 2 ) (tree dia- 
gram) and at 0(Q A ) (loop diagrams): 
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with k = k 2 = —k\. The one-pion-exchange contribution 
including the vertex renormalization (loop corrections to 
the vertices) reads 



Wtf = - 



(B.5) 



where t = fcp — k 2 ~ — k 2 , and F± (t) is the isovector 
Dirac form factor of the nucleon electromagnetic cur- 
rent. The phenomenologically determined Fi(t) is of 
the dipole type 



FY {t) 
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A 2 -t 
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with A = 840 MeV. The HBxPT expression for F^ (t) 
up to one-loop accuracy is given by 



APPENDIX B: AXIAL CHARGE OPERATORS 

As stressed in the main text, the axial charge opera- 
tors are chiral-protected. Since the axial-charge opera- 
tors up to 0(Q 4 ) have already been described in detail 
in Ref. [^5| , we only briefly recapitulate what is directly 
relevant to the present work. The leading one-body A 
operator is kinematically suppressed because of the 75 
matrix. Correspondingly, in xPT, the A operator at 
order C(Q 4 ) appears as a 1/mjv term. The leading cor- 
rection to the one-body axial-charge operator comes from 
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The loop functions, if's, will be specified below. The 
constant is determined by the nucleon charge radius ,[^| 



16 Note that this has nothing to do with the c's that 
appear in the GT operators. This confusing situation arises 
due to the lack of a unified system for labelling the coefficients 
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We should mention here that M.\-k in Eq^BJj) contains 
both <D(Q 2 ) and 0(Q A ) contributions. The 0{Q 2 ) con- 
tributions can be obtained by replacing (t) by 1, while 
(F^(t) — 1) is responsible for the 0(Q 4 ) contributions. 
The two-pion-exchange contributions, which are also of 
0(Q 4 ), are given by : 
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where and are unknown parameters^. The total 
two-body axial-charge operator is the sum of Eqs. (B.5) 
and (HI): 



= w£P + w£P (b.io) 

The loop functions K's in the above are defined as 
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The integrations over z can be done analytically, result- 
ing in 
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with 
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APPENDIX C: REGULARIZATION AND THE 
CUTOFF 

1. Fourier transform 

Since SNPA employs coordinate-space representation, 
we need to Fourier transform the momentum-space ex- 
pressions. In doing so, we must impose a cutoff to reg- 
ularize the integral. The cutoff introduced here typi- 
cally represents a scale that divides the low-energy de- 
grees of freedom (which we choose to include explicitly) 
and the high-energy degrees of freedom (which we in- 
tegrate out). How to implement cutoff into the the- 
ory is not unique, but physics should be independent 
of methods used insofar as the calculation is done con- 
sistently. This is a statement of renormalization group 
invariance. We now describe a particular cutoff scheme 
to be used here.|^ For the n^-body current in momen- 
tum space, A" 2 ... n = A" 2 ...„(fci, fc 2 , • • • , fe n ), define its 
"Fourier transform" as 
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where 5a (k 2 ) is a regulator with a cutoff A, and the 
factor (27r) 3 <j( 3 ) (q + k\ + fc 2 + ■ ■ ■ + k n ) comes from mo- 
mentum conservation. We employ here a regulator of the 
Gaussian type p| 



Sa (k 2 



exp 



2A 2 J 
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For a on e-body operator, the regulator plays no role, see 
Eq.flAj). Now for the two-body current, Eq.(|cl|) gives^ 



in the chiral Lagrangian. In a chiral expansion £ — £o + 
Ci + £i + • • ■, it was once common to use the letter b (c) to 
label the coefficients appearing in £\ (£2)- Later, however, 
it has become more common to use the letter c to designate 
the coefficients that feature in Ci, resulting in unfortunate 
confusions. In short, the above c§ is a coefficient in £2, while 
the c's appearing in the expression for the GT operator belong 
to £\. 

17 In fact, these are the only unknown parameters up to 
N 4 LO. It turns out that due to the chiral filter mechanism, 
one can drop these terms when working to N 4 LO for the GT 
operator. 



We illustrate our method for A, but the same method is 
used for the other currents as well. 

19 Using a Gaussian cutoff can in principle upset the chiral 
counting even when graphs up to a definite chiral order (say v) 
are considered. This counting mismatch, however, occurs at 
an order higher than p, and furthermore the error committed 
is likely to be minimized by the regularization scheme we have 
adopted for the d counter term. 

20 From here on, we shall specialize ourselves to the threshold 
limit, q = 0. 
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To simplify subsequent expressions, we define the follow- 
ing functions: 



Voa (m,r) 
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These functions become the ordinary delta and Yukawa 
functions when A goes to infinity. We also use the ab- 
breviation yJ A (r) = yoA( m u-, 0, an d similarly for Ui A (r) 
and ( r ) ■ The regularized delta and Yukawa functions 
read 



5< \\ r ) = 77~TX exp 

(47Tj 2 
1 m 2 1 

— (r <-* — r)] 



AV 
4 

- mr erfc 



Ar m 
T+A 



(C.5) 



(C.6) 



We are now ready to write down the i/=3 two-body 
axial-vector current, Eq.(A.£) in coordinate space: 



(fi x T 2 ) a (p\ cr 2 • fi2 +p>i • fi 2 ). (CIO) 



Note that Op a is completely determined by Lorentz in- 
variance. In terms of these seven operators, we can write 
all the two-body currents (including v—i and v=A con- 
tributions) as 

^12=- E [^(^i 2 )^ Q +^(ri 2 )^ Q 

0=+,-,X 

2 



2m N f 2 



2m N f 2 



E 



L0=+,-,x 



(C.ll) 



We have separated out here the part proportional to 
S^\r). The dimensionless parameters d,Q are given by 
d\ ji and (higher order) loop contributions as 



d+ = d x + -c 3 
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g A m N m m 
32tt 
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A 12 (1tt) = 
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where the superscript (i, j) are particle indices, ri 2 = 
|t*i2 I and ri 2 = In the above equations, we have 

defined the following two-body spin-isospin operators 



V 1 cr "■' 







(C.9) 



where = (+,—, x) and 



The functions Fq' T in Eq.( C.ll ) are given by 
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where fc = |fe|, Di = Di(k) and 
■d 3 k 



{/(fc 2 )} FT (r)^ j^Sl(k 2 )e-^f(e), 

{/(^ 2 )}pt U = 'Ul*J$F sl{k2) e " kr m - 

The explicit results of Fourier transformation of Di{k) 
and k 2 Di(k) are given by ("—>•" denotes Fourier transfor- 
mation) : 
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Next we turn to the three-body currents of Eq. ( A.lTj ). 
The Fourier-transformed three-body current has the form 
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with 
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The calculation of JJ 23 is rather involved. We may start 
with exploiting the identity 



(2^) 3 ^ 3 )(fc 1 +fc 2 + fe 3 ) = J d 
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to arrive at 
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1 123 — 
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where 



A= V2A. 



This representation is nicely transparent, and the result- 
ing integrand is non-oscillatory and rapidly damping. 
The remaining integration can be done by means of a 
Monte Carlo simulation with Metropolis random walks. 



2. The parameter d R 

Up to N 3 L0, unknown parameters occur only in A. 
At N LO, several unknown parameters appear in both 
V and Aq, but no new parameters appear in A. By the 
chiral filter argument, one can ignore the N LO terms 
in both V and Aq while going to N 4 LO in A. Thus, 
up to N 4 LO, the only genuinely unknown parameters re- 
side in A. The crucial observation is that up to N A LO, 
there is effectively only one constant d R that governs the 
GT amplitudes of all the cases under consideration. The 
argument goes as follows. 



The two parameters, d\ and d 2 , and the four A's can be 
combined into three unknown parameters, d± tX , that re- 
flect short-range physics. It is the Fermi-Dirac statistics 
that reduces the number of unknowns from three to one. 
To see this, let S CT and S T be the exchange operators in 
spin and isospin spaces, respectively; S CT = ^(1+ax -cr^), 
and S T = |(1 + ti ■ t 2 ). An explicit calculation gives 
the identity tj\ x <x 2 = i{<j\ — cr 2 )'E J a : and likewise for 
t\ x T2- Now, the Fermi-Dirac statistics requires that 
S r S°'S r = — 1, where H r is the Majorana exchange op- 
erator that exchanges the orbital coordinates, r± and r 2 . 
As a result, 



(C.20) 



When multiplied by the delta function S^ 3 \r), the oper- 
ators are non-vanishing only for the L — states, which 
then implies S + T = 1. Acting on L = states, O l l a 
is identically zero, since either spin or isospin must be 
equal to zero. Furthermore, the L = states are eigen- 
states of the operator S r with eigenvalue 1, so that 0!l a 
becomes identical to O'/ . Thus we are left with only one 
unknown parameter, d = d- + d x . 

The above argument is not strictly valid for the cut- 
off delta function S A (r), which has a finite (albeit very 
small) range, ~ A -1 . However, deviations from the ordi- 
nary delta function case is higher order in chiral counting 
and hence can be ignored. 



APPENDIX D: CALCULATION UP TO N^LO 

We have derived all the weak currents up to N 4 LO. As 
Table | indicates, loop contributions start at N 4 LO. Loop 
corrections in the vector currents (both V and V°) can 
be safely ignored, since even their leading single-particle 
terms are suppressed relative to the axial current. It 
turns out that the loop diagrams in A are all finite and 
hence need no regularization although there are finite 
counter terms that should be taken into account. On the 
other hand, the loop diagrams in Aq do have divergences 
and need to be regularized. To derive the momentum 
space expressions for the currents given above, we have 
employed the dimensional regularization. This is not 
quite congruous with the cutoff regularization adopted in 
going from momentum to coordinate space. Meanwhile, 
using a cutoff regularization in evaluating loop graphs is 
a delicate matter, since that might endanger chiral sym- 
metry; with the use of a cutoff regularization one might 
need chiral-symmetry-breaking counter terms in order to 
satisfy the Ward identities. We have not yet investigated 
whether the dimensional regularization as used here pre- 
serves chiral symmetry, and hence we cannot say at this 
point whether our coordinate space operators at N 4 LO 
arc fully consistent. However, this problem does not arise 
if we limit ourselves to N 3 LO, for up to this order there 



are no loop contributions. The relevant 2-body currents 
in coordinate space are: 



Vi2(r) 

A\ 2 {r) 
A 12 (r) 



iiA.it x 



- 9 A 

|(o + 

9A 

2m N flr 



T x --0 



—k ■ f <r_ 

2 



O-) 

+ r 



c 4 + - 



O 



-T x (p x <T 2 



2/2AM 

r +p 2 cri ■ r)yl K (r) 



4 3) (r)^0> 



(D.l) 



with 



d 3 k 



■ e a- 1 e 



1. 



k 2 

d 1 



vUr) = -r&yfoir) and yJ A (r) = ^(r), 



where Oq 



1 = r m cr k 



n-k = 

0"0' ^0 — 



0% and 



d x + 2d 2 



7C3 



7C4 



(D.2) 



The derivative operators, p t (I — 1, 2) in eq. (f22|) , should 
be understood to act only on the wave functions. 

Since dn accompanies the regularized ^-function 
<5 A (r), its contribution depends on A rather strongly. 
However, the renormalization-group invariance of EFT 
requires that this sensitivity to A should be compensated 
by the contributions of the remaining terms. Since the 
single-particle piece of A has no A dependence, and since 
all the currents other than A have only weak A depen- 
dence, this compensation must occur between the finite- 
range two-body GT and the regularized delta-function 
term. This has been indeed verified in our calculation 
over a wide range of A (500 MeV - 800 MeV) although, 
as mentioned above, the results for the 800 MeV cutoff 
should be viewed with caution. 
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